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Executive  Summary 

The  U.S.  Navy  is  currently  exploring  the  adoption  of  novel  ship  design  concepts  that  support  off¬ 
shore  basing  and  littoral  combat  operations:  examples  include  the  Joint  High-Speed  Vehicle 
(JHSV)  and  Littoral  Combat  Ship  (LCS)  programs.  As  part  of  the  JHSV  and  LCS  programs, 
light-weight  materials  ( e.g .,  aluminum)  and  non-conventional  hull  forms  ( e.g .,  eatamaran-style 
hulls)  are  being  evaluated  for  adoption.  Compared  to  steel  mono-hull  ships,  the  U.S.  Navy  has 
comparatively  less  design  and  operational  experience  with  high-performance  aluminum  vessels. 
Hence,  a  dire  need  exists  for  the  development  of  comprehensive  life-cycle  monitoring  and 
assessment  strategies  that  can  track  the  performance,  state  of  health  and  expected  lifespan  of 
operational  high-speed,  aluminum  ships  exposed  to  demanding  seaway  loads.  In  response  to  this 
life-cycle  monitoring  framework,  the  Naval  Surface  Warfare  Ccnter-Carderock  Division 
(NSWCCD)  has  proposed  a  comprehensive  research  program  in  improving  the  characterization 
of  the  structural  materials  for  ship  reliability.  The  life-cycle  framework  includes  the 
development  of  effective  structural  health  monitoring  (SHM)  methodologies  to  support  early 
detection  of  structural  failures  (e.g.,  fatigue  failure). 

This  research  project  represents  a  preliminary  investigation  into  a  comprehensive  analytical, 
numerical  and  experimental  framework  for  the  monitoring  and  life-cycle  assessment  of  the 
structural  integrity  and  performance  of  aluminum  hull  structure.  Aluminum  plate  specimens  arc 
designed  and  fabricated  to  facilitate  the  investigation  of  system  identification  and  damage 
detection  methodologies,  both  of  which  are  key  components  of  future  life -cycle  analyses.  The 
design  of  the  plate  specimens  is  intended  to  include  the  geometric  complexity  commonly  found 
in  aluminum  ship  hull  structures.  A  powerful  Bayesian  probabilistic  model  updating  framework 
is  applied  to  validate  its  potential  applicability  for  damage  identification  around  critical  weld 
zones  where  fatigue  failure  is  likely  to  initiate.  Extensive  numerical  simulation  and  experimental 
testing  has  been  conducted  to  validate  the  model-based  approach  to  crack  damage  detection  in 
aluminum  plate  structures.  Future  efforts  include  continuing  collaboration  with  the  researchers 
at  NSWCCD  in  the  fatigue  testing  of  aluminum  ship  hull  components  and  continued 
development  of  the  model-based  damage  detection  methodology  as  an  integral  part  of  a  more 
general  lifc-cyclc  design,  monitoring  and  maintenance  framework.  The  long-term  benefit  of  the 
proposed  SHM  strategy  is  the  reduced  total-cost -of-ownersh ip  for  aluminum  littoral  combat 
ships  in  addition  to  more  resilient  ship  structures  for  use  in  demanding  seaway  environments. 
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1.  Objectives 

The  U.S.  Navy  is  currently  exploring  the  adoption  of  novel  ship  design  eoneepts  that  support  off¬ 
shore  basing  and  littoral  eombat  operations;  examples  inelude  the  Joint  High-Speed  Vehiele 
(JHSV)  and  Littoral  Combat  Ship  (LCS)  programs  (Hess  2007).  As  part  of  the  JHSV  and  LCS 
programs,  light-weight  materials  (e.g.,  aluminum)  and  non-eonventional  hull  forms  (eg., 
catamaran-style  hulls)  are  being  evaluated  for  adoption.  In  addition  to  the  material  properties 
being  high  corrosion  resistant,  the  use  of  aluminum  alloys  in  the  design  of  naval  structures  offers 
the  benefit  of  light-weight  ship  constructions  that  favor  high  speeds  at  sea.  However,  eompared 
to  steel  mono-hull  ships,  the  U.S.  Navy  has  comparatively  less  design  and  operational  experience 
with  high-performanee  aluminum  vessels.  Henee,  a  dire  need  exists  for  the  development  of 
comprehensive  life-cycle  monitoring  and  assessment  strategies  that  can  measure  the  performance, 
state  of  health  and  expeeted  lifespan  of  operational  high-speed,  aluminum  ships  exposed  to 
demanding  seaway  loads.  Current  practice  in  engineering  design  and  maintenance  are  primarily 
oriented  towards  serviceability  and  safety  but  typieally  exclude  consideration  of  possible  damage 
in  the  structure.  For  life-eyelc  design,  it  is  critically  important  to  go  beyond  the  basie  load 
scenarios,  standard  material  behavior,  and  serviee  and  safety  criteria.  Methods  for  lifetime- 
oriented  design  must  be  able  to  describe  damage  initiation,  damage  evolution,  and  remaining  life 
resistance,  as  well  as  to  estimate  the  reliability  and  performance  of  a  structural  system  ineluding 
those  with  damage  (Stangenberg  et  al.  2009). 

The  use  of  aluminum  poses  a  number  of  technical  challenges  for  the  naval  engineering 
community  including  a  higher  incidence  of  fatigue-related  cracks.  Fatigue-related  damage  in 
aluminum  alloys  often  appear  as  widespread  micro-cracks.  High-speed  aluminum  hulls  can 
remain  in  operation  even  after  the  initiation  of  micro-eraeks  bceause  of  the  duetile  mechanical 
characteristic  of  the  material.  However,  fatigue  cracks  in  aluminum  can  grow  quiekly  resulting 
in  rapid  deterioration  of  the  hull.  Therefore,  one  key  component  to  the  lifc-cyclc  monitoring 
framework  is  the  early  detection  of  fatigue-indueed  eraeks.  Early  detection  of  potential  damage 
to  the  aluminum  ship  structure  not  only  enhances  the  effectiveness  of  maintenance  strategies, 
but  also  prevents  the  catastrophic  failure  of  the  hull.  Furthermore,  monitoring  the  integrity  of  the 
aluminum  hull  can  provide  valuable  information  that  ean  be  used  to  aeeurately  estimate  the 
residual  life  of  hull  components. 

The  Naval  Surface  Warfare  Center  -  Carderock  Division  (NSWCCD)  has  proposed  a 
comprehensive  research  program  in  improving  the  characterization  of  structural  aluminum  alloy 
materials  for  ship  reliability.  The  NSWCCD  effort  places  a  particular  emphasis  on  fatigue 
cracking  in  aluminum  welded  connections  and  assemblies  common  to  LCS  ship  designs.  In 
collaboration  w'ith  researchers  at  NSWCCD,  this  study  represents  a  preliminary  investigation  of 
a  framework  for  the  monitoring  and  performance  assessment  of  aluminum  hull  components. 
Specifically,  aluminum  plate  specimens  have  been  designed  and  built  w  ith  welded  assemblies  to 
facilitate  the  investigation  of  a  system  identification  and  damage  detection  methodology.  This 
report  describes  a  model-based  damage  detection  methodology  for  autonomous  structural  health 
monitoring  of  aluminum  hull  structures  monitored  with  sensors.  The  area  of  damage  in  the 
structure  is  estimated  by  comparing  the  structural  characteristics  of  the  “true”  structure  (damaged 
or  undamaged)  and  hypothesized  (or  “trial”)  finite  element  method  (FEM)  models.  The 
probability  associated  with  a  hypothesized  damage  state  (e.g.,  location  and  size)  is  evaluated 
through  the  calculation  of  an  error  between  the  “true”  and  “trial”  models.  The  probable  damage 
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area  (i.e.,  fatigue  eraek  path)  is  identified  by  repeatedly  applying  a  Bayesian  inference  algorithm 
(Sohn  and  Law  2002).  This  inference  framework  has  the  advantage  of  providing  a  probabilistic 
basis  for  guiding  the  selection  of  trial  models  with  the  aim  of  substantially  reducing  the 
computational  demands  associated  with  exhaustive  searches.  The  Bayesian  probabilistic  model¬ 
updating  framework  is  explored  in  this  study  for  its  potential  applicability  to  damage 
identification  around  critical  weld  zones  on  standard  aluminum  plate-stiffener  specimens. 
Specifically,  extensive  numerical  simulations  and  experimental  tests  have  been  conducted  to 
examine  the  model-based  Bayesian  approach  for  the  structural  health  monitoring  of  aluminum 
hull  components  that  are  vulnerable  to  fatigue  failure. 


2.  Theoretical  Background  and  Algorithm 
2.1.  Model  Updating  for  Damage  Detection 

A  general  class  of  damage  detection  algorithms  is  based  on  the  concept  of  model  updating 
(Doebling  et  al.  1 998).  In  model  updating,  parameters  in  an  analytical  model  (typically  based  on 
physical  principles)  arc  varied  until  the  model  approximates  the  behavior  of  the  true  observed 
system.  To  identify  an  optimal  model,  the  outputs  from  the  observed  system  (namely,  sensor 
measurements)  are  used  to  evaluate  a  pre-defined  objective  (or  error)  function.  An  appropriate 
objective  function  is  one  that  takes  into  aeeount  the  fundamental  behavior  of  the  system  in  both 
its  damaged  and  undamaged  states,  yet  is  compatible  with  the  experimental  data  available. 
Changes  in  the  model  parameters  are  then  correlated  to  the  (damaged  versus  undamaged) 
condition  of  the  structure. 


Consider  a  simple  plate  structure.  The  governing  equation  describing  the  dynamic  behavior  of  a 
vibrating  plate  can  be  written  in  the  following  form  (Chakravcrty  2008): 


DV4w+  ph 


q(x,y,t) 


Eq.  1 


where  D  is  the  plate  flexural  rigidity  and  is  defined  as: 


D  = 


Bf_ 

12(1  -v2) 


Eq.  2 


Here,  w  denotes  the  vertical  displacement  of  the  plate,  q(x,y,t)  is  the  normal  load  distribution 
function  on  the  top  of  the  plate,  p  is  the  plate  density  (mass  per  unit  area),  h  is  the  plate 
thickness,  E  is  the  Young’s  (clastie)  modulus  and  v  is  the  Poisson  ratio  of  the  plate  material. 
The  in-plane  coordinates  of  the  plate  are  denoted  as  x  and  y  while  t  denotes  time.  To  model  the 
plate  under  eomplex  loading  conditions  or  when  the  plate  possesses  non-uniform  material 
properties,  the  finite  element  method  (FEM)  ean  be  used  as  an  effective  analytical  method 
(Hughes  2000). 

For  a  finite  element  model  with  N  elements,  the  elastic  moduli  for  the  elements  ean  be  denoted 
as: 
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E  [  E] ,  E2 ,  •  •  • ,  E, ,  Em  Ej+n ,  •  •  • ,  En_ , ,  En  J 


Eq.  3 

For  modeling  purposes,  this  project  models  fatigue  crack  damage  as  a  change  in  the  flexural 
rigidity  or  stiffness  of  the  plate,  which,  in  turn,  is  reflected  by  a  reduction  in  the  clastic  modulus 
of  the  damaged  elements.  If  the  model  contains  n  damaged  elements,  their  clastic  moduli 
{^l+l  ’  ‘  ‘  ‘  ’  Ei+ Ei+n  ^  are  replaced  by 

Ed={K  A+v—’ki+jEi+j>—’k,+nE,+n}  Eti- 4 

where  kUj  <1  (j=l,...,n)  denotes  the  reduction  factor  for  the  (i+j)'h  damaged  element. 

Model  updating  as  the  primary  means  of  performing  damage  detection  is  essentially  a 
combinatorial  optimization  (CO)  problem  that  seeks  to  Find  the  optimal  set  of  elastic  moduli  that 
minimizes  a  defined  objective  function  by  comparing  the  FEM  model  output  and  the  actual 
system  output  derived  from  sensor  measurements.  Assuming  an  linear  variation  in  kUj  .the  CO 

problem  is  defined  by  an  infinitely  large  state  space  that  is  impossible  to  exhaustively  search  to 
identify  the  best  model  parameters.  In  addition,  the  objective  functions  typically  used  in  CO  arc 
rarely  eonvex  function;  rather,  they  are  defined  by  a  family  of  local  minima  one  of  which  is  the 
true  global  minimum  (i.e.,  the  optimal  model).  As  a  result,  gradient  decent  methods  arc  prone  to 
converge  to  local  minima  which  can  be  far  from  the  true  global  minimum.  As  a  result,  CO 
problems  are  extremely  complex  problems  to  solve.  In  this  study,  a  probabilistic  Bayesian 
methodology  is  adopted  to  balance  computation  effort  with  the  accuracy  of  the  model  updating 
solution. 

2.2.  A  Flexibility-Based  Objective  Function 

Conceptually,  damage  detection  is  to  identify  changes  (i.e.,  damage)  in  a  structure  by  using 
specific  structural  characteristics  as  the  basis  for  evaluation.  Vibration  data  are  among  the  most 
common  sensor  information  easily  obtained  from  dynamic  excitations.  Vibration-based  damage 
identification  methods  using  mode-shape  information  have  been  proposed  (Docbling  et  al.  1998). 
One  approach  is  to  detect  damage  directly  based  on  the  changes  in  structural  eharactcristics 
between  damaged  and  undamaged  structures.  Examples  of  structural  characteristics  include 
curvature  mode  shape  (Pandcy  et  al.  1991),  flexibility  parameters  (Pandcy  and  Biswas  1995, 
Bernal  2006),  strain  and  modal  strain  energy  (Stubbs  et  al.  1995;  Zonta  and  Bernal  2006). 
Although  current  research  in  damage  detection  has  made  substantial  progress,  the  methods 
developed  so  far  are  primarily  restricted  to  simple  beam  or  frame  structures  with  a  limited 
number  of  degrees-of-freedom.  In  contrast,  continuum  systems  such  as  plate-like  structures  pose 
many  challenges  for  current  methods  because  of  the  large  number  of  degrees-of-frecdom  that  arc 
required  to  properly  model  the  structure.  Additionally,  it  is  not  realistic  to  place  a  sufficient 
number  of  sensors  to  measure  the  system  response  at  all  of  its  dcgrecs-of-frccdom.  An 
alternative  approach  is  to  construct  a  system  model  (such  as  a  finite  element  model)  of  a  targeted 
structure  using  measurement  data.  This  “model-based”  approach  updates  the  system  model  by 
modifying  the  structural  properties  of  the  elements  until  the  model  resembles  the  dynamic 
characteristics  estimated  from  sensor  information.  The  objective  is  thus  to  define  an  objective 
(error)  function  and  to  compute  the  difference  (error)  between  the  trial  finite  clement  model  and 
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the  target  structure.  Damage,  if  any,  is  identified  and  revealed  through  subtle  changes  in  the 
model  parameters  extracted  during  the  model  updating  process. 

For  the  plate  problem,  our  study  indicates  that  the  use  of  flexibility  matrices  constructed  based 
on  modal  properties  ( i.e .,  modal  frequencies  and  mode  shapes),  as  opposed  to  the  direct  use  of 
modal  properties,  is  one  appropriate  choice  for  the  objective  function  (Kurata  et  al.  2010).  The 
inverse  relationship  between  the  flexibility  matrix  and  the  square  of  the  modal  frequency  renders 
the  flexibility  matrix  less  sensitive  to  high  frequency  modes  which  are  difficult  to  identify  in 
monitoring  data  (due  to  the  electrical  noise  inherent  to  the  data  collection  process).  This  unique 
characteristic  allows  for  the  inclusion  of  lower-order  modes  in  a  truncated  flexibility  matrix. 
This  feature  has  attracted  many  researchers  to  explore  flexibility  as  a  core  element  in  developing 
structural  damage  detection  algorithms  (Bemel  2006;  Gao  and  Spencer  2002;  Pandcy  and  Biswas 
1994;  Pandey  and  Biswas  1995;  Zonta  and  Bernal  2006).  In  the  so-called  “flexibility-based 
approach”  (FBA)  described  herein,  the  objective  (error)  function  is  expressed  in  terms  of  the 
difference  between  the  flexibility  matrices  that  correspond  to  the  true  system! measured)  and  trial 
models  (finite  element). 


When  the  mode  shapes  are  mass-normalized  (i.e.,  <p'  Mfa  =  /),  the  flexibility  matrix  F  can  be 
expressed  in  terms  of  the  system  modal  properties  as  follows; 


F  =  Eq.5 

1-1  C 0 ; 

where  M  is  the  mass  matrix,  fa  -\j\fa2 . . ,faN~^  is  the  mode  shape  matrix,  Q  =  diag(ar)  is  the 
spectral  matrix  consisting  of  the  square  of  the  modal  frequencies,  coi ,  and  N  is  the  number  of 
degrees-of-freedom  in  the  system.  The  mass-normalized  modal  vector,  faj  ,  is  related  to  the 
arbitrarily  sealed  mode  shape,  fa.  as: 

fa  =  fadi  Eq.  6 

where  d.  =  is  a  mass  normalization  constant  for  the  i'h  mode. 

# AM 

Suppose  only  a  few  (typically  lower)  modes  are  available,  then  a  truncated  flexibility  matrix  is 
obtained  as: 


=£« 


'  d. v 


j 


where  n  denotes  the  number  of  modes  available. 


Eq.  7 
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Let  us  define  the  difference  ( A Ftrun )  between  the  flexibility  matrices  of  the  true  (damaged) 
structure  and  the  trial  FEM  model  as: 


A  f  =  p,n,e  _  f?'rial 
trim  *  trun  1  trun 


Eq.  8 


When  a  trial  FEM  model  reasonably  resembles  the  damaged  structure,  the  difference  in 
flexibility  matrices  is  close  to  zero  (exactly  zero  if  there  is  no  measurement  noise  nor  modeling 
error).  The  scalar  magnitude  on  the  difference  in  the  flexibility  matrices  can  be  measured  by 
calculating  the  Frobenius  norm  of  A Ftnm : 

K-L-i/LI,**  E<i-9 

which  vanishes  when  all  matrix  elements  x,y  in  A are  zero  (meaning  the  FEM  model  perfectly 

matches  the  observed  structure).  The  difference  in  the  flexibility  matrices  can  also  be  further 
decomposed  into  singular  values  by  singular  value  decomposition  (SVD): 

A F,nm=USVT  Eq.  10 


where  VT  and  U  arc  matrices  of  singular  vectors,  and  S  is  the  diagonal  matrix  whose  elements 
are  the  singular  values,  st .  Since  the  Frobenius  norm  is  invariant  under  unitary  transformation, 
the  SVD  of  the  A yields 


+  .?2  +  •••+$# 


Eq.  11 


where  R  is  the  rank  of  A .  Thus,  the  objective  is  to  search  for  a  trial  FEM  model  that 
minimizes  the  Frobenius  norm  shown  in  Eq.  1 1 . 

Since  the  mode  shapes  experimentally  obtained  arc  arbitrarily  scaled,  the  mass  normalization 
constants  (c/.)  are  required  a  priori  to  properly  compute  the  flexibility  matrix  of  the  real,  true 

structure.  One  approach  to  extracting  the  mass  normalization  constants  is  based  on  testing  the 
structure  with  a  perturbed  mass  matrix  (by  adding  a  known  mass  at  a  certain  location)  and 
examining  the  sensitivity  of  the  eigenvalues  (Brinker  and  Anderson  2002;  Parloo  et  al.  2002).  In 
this  study,  mass  normalization  constants  estimated  using  the  FEM  model  are  updated  at  each  trial 
and  applied  to  the  experimentally  obtained  (and  not  mass  normalized)  mode  shapes.  That  is,  the 
constants  are  extracted  by  comparing  the  mass-normalized  mode  shapes  and  displacement- 
normalized  mode  shapes,  where  both  mode  shapes  are  available  as  options  in  many  commercial 
FEM  programs  ( e.g .,  ABAQUS). 
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2.3.  Bayesian  Probabilistic  Approach 

The  model  updating  procedure  adopted  in  this  study  is  based  on  a  Bayesian  probabilistic 
approach  which  utilizes  the  parameters  measured  or  estimated  from  a  series  of  collected 
vibration  signals  or  data.  Unlike  a  deterministic  optimization  formulation,  the  state  spaee  search 
must  reflect  the  relative  degree  of  belief  on  the  estimates  of  the  optimal  subset  (i.e.,  Ed  in  this 

ease).  Let  M  denote  the  hypothesized  damage  states  of  the  model.  The  calculation  of  the  error 
that  exists  between  the  “true”  structure  and  the  FEM  “trial”  model  is  based  on  the  measured  or 
estimated  structural  parameters.  The  updated  estimates  on  the  damage  of  the  structure  are 
expressed  as  the  posterior  distribution  based  on  Baye’s  rule  as  follows: 

p(M  |  s)  =  oc  p(s  |  M)p(M)  Eq.  12 


where  p(M  \  s)  is  the  posterior  distribution  function  for  a  hypothesis  M  given  the  measured  or 
estimated  parameters,  s.  p(s  \  M)  is  termed  the  likelihood  function,  p(M)  is  the  prior 
probability  of  the  hypothesis  and  p(s)  is  treated  here  as  some  normalization  constant.  By 
collecting  the  likelihood  function,  the  posterior  distribution  p(M\s)  would  progressively 
become  a  better  estimate  than  the  prior  probability  p(M)  as  the  process  goes  on.  For  instance, 
if  the  initial  estimate  of  the  probable  area  of  damage  is  assumed  to  be  uniformly  distributed  over 
the  structure  before  the  model  updating  procedure,  the  most  likely  damaged  areas  are  revealed 
with  relatively  higher  posterior  estimates  by  the  repeated  application  of  the  Bayesian  inference 
algorithm. 

The  selection  of  the  most  probable  events  from  all  conceivable  possibilities  using  the  Bayesian 
probabilistic  approach  ean  be  systematically  implemented  using  various  “optimal”  search 
methods;  otherwise  a  random  search  in  the  state  spaee  beeomes  a  computationally  intraetable 
task.  To  sample  the  posterior  parameter  distribution,  the  Bayesian  approaches  implemented  with 
the  Markov  Chain  Monte  Carlo  (MCMC)  method  and  genetie  algorithms  (GA)  have  been 
reported  with  successful  performance  in  detecting  structural  damage  detection  (Vanik  et  al. 
2000;  Cheung  and  Beck  2009;  Stull  et  al.  2009;  Niehols  et  al.  2009).  However,  these  approaches 
are  computationally  expensive  and  they  do  not  guarantee  convergence  to  an  optimal  solution.  To 
reduce  the  computational  effort,  the  Bayesian  damage  detection  algorithm  proposed  herein  is 
enhanced  with  a  braneh-and-bound  (BB)  search  technique  where  the  state  spaee  searched  is 
systematically  narrowed  through  enumeration  and  pruning  of  candidate  model  solutions  (Sohn 
and  Law  2002). 

The  BB  technique  is  a  general  search  method  originally  developed  for  discrete  optimization 
problems  and  is  a  powerful  technique  for  controlling  the  size  of  a  search  spaee  used  in  model 
updating  (Norkin  et  al.  1998).  As  illustrated  in  Figure  1,  the  BB  algorithm  initially  starts  its 
search  from  some  random  subspaee  (i.e.,  leaf  nodes  associated  with  Branch  1  in  Figure  1(a)).  At 
each  step,  the  algorithm  takes  an  additional  sample  (i.e.,  an  additional  element  in  the 
hypothesized  subset)  at  each  leaf  node  and  aims  to  improve  its  estimates  (this  process  is  called 
branching).  As  the  search  proceeds,  the  branches  associated  with  large  “errors”  are  eliminated 
from  further  consideration  and  the  search  is  bounded  by  evaluating  the  remaining  branches  (this 
process  is  called  pruning).  Figure  1(b)  schematically  shows  the  application  of  the  BB  technique 
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Branch  1 


(a)  Concept  of  the  branch-and-bound  search  scheme 

"Current  estimate"  "Next  candidates" 

Crack  damage 

^rFP  rFF*  rFF  ‘  rffl  FEmodel 

(b)  Candidates  for  probable  damage  state 
Figure  1.  Application  of  the  branch  and  bound  technique  for  model-updating 


Estimated  location 
of  crack  initiation  . 


for  FEM  updating  of  a  plate-type  strueture.  Here,  the  eraek  is  assumed  to  damage  the  entire 
section  of  the  plate  in  the  thickness  direction  and  represented  in  the  model  by  a  set  of  finite 
elements  with  reduced  clastic  modulus  for  simplicity.  The  collection  of  survived  leaf  nodes  (or 
“trial”  model  with  small  error)  is  mapped  into  a  damage  map  at  every  branching  process. 
Therefore,  the  probable  damage  area  in  a  system  can  be  systematically  updated  by  implementing 
the  BB  technique  in  the  Bayesian  formulation.  The  computational  effort  and  accuracy  of  the 
estimates  for  the  optimal  subset  ( i.e .  Ed )  can  be  controlled  by  selecting  an  optimal  pruning  rate 
at  each  branch. 

2.4.  Summary  of  the  Model  Based  Bayesian  Damage  Detection  Approach 

The  model-based  structural  health  monitoring  framework  consists  of  three  phases  as  shown  in 
Figure  2.  The  first  phase  is  the  construction  of  the  FEM  model  with  structural  properties  and 
boundary  conditions  close  to  that  of  the  target  structure.  In  the  second  phase,  the  monitoring 
system  acquires  information  about  the  status  of  the  target  structure  through  a  network  of  sensors 
installed  in  the  structure.  Onee  the  variation  in  the  modal  properties  exceeds  a  pre-defined 
threshold  value,  the  third  phase  for  damage  detection  begins.  The  model-based  damage  detection 
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algorithm  computes  the  error  associated  with  an  initial  set  of  hypotheses  using  an  objective 
function.  Engineering  judgment  based  on  knowledge  about  common  damage  patterns  can  be 
very  useful  to  account  for  the  selection  of  the  initial  set  of  hypotheses;  for  instance,  heat  affected 
zones  (HAZ)  around  the  hull  welds  or  notches  along  component  edges  are  good  initial  candidates 
for  the  hypothetical  damage  initiation  areas.  The  branching  process  adds  one  more  damaged 
element  to  each  hypothesis  following  the  branching  rule  shown  in  Figure  1(b).  The  hypotheses 
with  relatively  large  error  arc  pruned  before  the  next  branching  process.  The  branching,  error 
computation  and  pruning  processes  make  a  finite  loop  until  the  damage  identification  algorithm 
terminates.  At  the  end  of  each  loop,  the  elements  surviving  in  the  current  and  previous  pruning 
processes  are  classified  and  stored  in  a  bin  based  on  their  error.  The  histogram-like  plot  of 
survived  elements  serves  as  a  damage  map  showing  the  probable  arca(s)  of  damage.  The  search 
proeess  is  terminated  when  the  deviation  between  the  current  and  last  damage  areas  converges 
within  a  predefined  tolerance. 
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Initial  model  setup  before 
damage  initiates 


Momtonng  the  structure  for 
triggering  event 


Crack  damage  detection 
algorithm 


Finite  element  model  of 
monitored  structure 


Update  finite  element  model 


Modal  properties  of 
monitored  structure 


Compare  modal  properties  of  FE 
model  and  monitored  structure 


Branching  process: 

Systematically  add  one  cracked  element 
to  candidate  FE  models 


Compute  objective  functions 
for  candidate  FE  models 


Pruning  process: 

Sort  candidate  FE  models  by  objective 
functions  and  prune  models  with  small  error 


Modal  properties  of 
monitored  structure 


Modal  properties  of 
monitored  structure 


Map  FE  models  with  small 
error  in  histogram  plot 


Compare  histograms  from 
current  and  prior  loops 


Histogram  plot 
from  prior  loop 


Figure  2.  Flowchart  of  the  model-based  Bayesian  damage  detection  system 
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3.  Numerical  Simulation  and  Experimental  Validation 

The  Bayesian  model-updating  procedure  is  applied  to  the  problem  of  damage  detection  on  a 
stiffened  aluminum  plate,  as  shown  in  Figure  3(a).  The  design  of  the  structural  plate  is  intended 
to  include  the  geometric  complexity  that  is  commonly  found  in  aluminum  ship  hull  structures. 
The  aluminum  plate  includes  an  area  with  high  stress  concentrations  due  to  the  presence  of  a 
welded  stiffener.  Such  areas  arc  the  likely  locations  for  fatigue-related  damage  initiation  (i.e., 
fatigue  crack).  Knowledge  of  this  fact  allows  one  to  customize  the  model  updating  algorithm  to 
prioritize  the  search  of  this  area  for  damage.  The  aluminum  plate  is  24-in  by  48-in  and  is  14-in 
thick.  In  addition,  the  plate  has  a  2-in  by  18-in  by  14-in  stiffener  plate  welded  to  it  off-center. 
The  base  plate  and  the  stiffener  plate  are  rigidly  welded  by  a  tungsten  inert  gas  (TIG)  weld. 
Three  different  crack  paths  initiating  from  the  HAZ  around  the  welded  stiffener  plate  are 
considered  as  example  damage  scenarios  as  shown  in  Figure  3(b)  and  denoted  as  crack  paths  1, 2 
and  3. 


2*  by  18",  1/4'  thick 

(a)  Schematics  of  the  test-bed  plate  structure 


(b)  Crack  paths  considered 
element 

sensorlocation^  20. ..23 


;.Nc 

H 

;  element  J 

:  *  element 

•  :  i _ 19 

:  "24.. .42 

* 

*X  1 

^  stiffener 


elements  probable  for/ 
damage  initiation 
(c)  Sensor  locations  and  probable  damage  zone 


Figure  3.  A  Test-bed  Plate  Structure 
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The  plate  structure  is  modeled  using  ABAQUS,  a  general-purpose  FEM  analysis  program.  The 
base  plate  and  stiffener  plate  are  modeled  with  4-node  reduced  integration,  doubly  curved  shell 
elements  with  hourglass  control  (S4RSIV)  (ABAQUS  2010).  These  plates  are  assumed  to  be 
rigidly  connected.  The  mesh  size  of  the  plate  elements  are  1  in  x  1  in.  The  clastic  modulus  E  of 
the  undamaged  elements  is  10300  ksi  while  for  the  cracked  elements,  the  clastic  modulus  is 
reduced  10f'  times  from  the  original  modulus  (Ed  =  Ex  10  *).  The  mass  density  of  the  aluminum 

is  assumed  to  be  2.489x10'  slug/in3.  The  Lanczos  frequency  analysis  method  in  ABAQUS  is 
employed  to  compute  the  modal  properties  ( i.e .,  modal  frequency  and  mode  shape)  of  the  base 
structure  as  well  as  the  structural  model  with  different  damage  scenarios  introduced.  Only  the 
first  five  modes  are  considered  in  the  numerical,  as  well  as  experimental  simulations. 
Additionally,  sensors  are  placed  (or  assumed  to  be  placed)  at  strategic  locations  as  shown  in 
Figure  3(c)  by  the  “x”  marks. 

3.1  Numerical  Simulations 

For  the  numerical  simulation  study,  the  test-bed  plate  structure  as  shown  in  Figure  3(a)  is 
employed  for  evaluating  the  model-based  damage  detection  algorithm.  Three  separate  cases  arc 
considered.  First,  damage  detection  of  the  plate  with  individual  cracks  and  initial  probable 
damage  zones  as  shown  in  Figure  3(b)  and  3(c),  respectively,  is  investigated.  Second,  similar  to 
the  first  damage  detection  case  with  the  three  individual  crack  paths,  the  plate  is  assumed  to  have 
another  welded  stiffener  and  the  initial  probable  damage  zone  is  expanded  to  include  other 
possible  damage  initiation  areas.  For  the  third  case,  damage  detection  with  multiple 
simultaneous  cracks  is  studied. 

Damage  Detection  of  Individual  Crack  Paths  with  Q.  Single  Stiffener:  The  modal  properties  for 
the  undamaged  plate  structure  and  the  damaged  plate  with  cracks  arc  obtained  first.  Table  1 
summarizes  the  modal  frequencies  obtained  for  the  undamaged  and  damaged  plates.  As  an 
example,  Figure  4  shows  the  mode  shapes  of  the  model  with  damage  introduced  as  crack  path  1. 
It  can  be  seen  that  the  first  and  third  modes  are  bending  modes,  the  second  and  fourth  are 
torsional  modes,  and  the  fifth  is  a  bending  mode  in  both  directions  of  the  plate.  Using  the 
computed  modal  frequencies  and  mode  shapes,  the  flexibility  matrices  of  the  damaged  models 
are  constructed  using  Eq.  7. 

Table.  1  Modal  frequencies  of  the  damaged  and  undamaged  plates 


Mode 

No  crack 
(Hz) 

Crack  1 

(Hz)  Diff  (%) 

Crack  2 

(Hz)  Diff  (%) 

Crack  3 

(Hz)  Diff  (%) 

1 

26.88 

26.52 

-1.34 

26.82 

-0.22 

26.86 

-0.05 

2 

41.76 

39.89 

-4.48 

41.60 

-0.39 

41.69 

-0.18 

3 

72.87 

70.89 

-2.72 

72.02 

-0.18 

72.82 

-0.08 

4 

94.93 

89.49 

-5.73 

94.88 

-0.05 

94.85 

-0.09 

5 

128.95 

126.93 

-1.57 

128.64 

-0.24 

128.93 

-0.02 
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Mode  1 


Mode  2 


Mode  3 


Figure  4.  Mode  shapes  of  the  baseline  structure  with  crack  I 

The  model-updating  procedure  (as  described  in  Figure  2)  is  employed  to  identify  the  individual 
eraeks  introduced  in  the  aluminum  hull  specimen.  Figure  5  illustrates  the  basie  proeess  in 
locating  eraek  path  1  shown  in  Figure  3(b).  Starting  from  probable  damaged  elements  in  the 
vieinity  of  the  weld  (i.e.,  HAZ  area  shown  in  Figure  3(e)),  at  eaeh  branching  step  the 
probabilistic  braneh-and-bound  seheme  proceeds  to  seareh  for  the  most  likely  set  of  damaged 
elements  based  on  the  objective  function  shown  in  Eq.  11,  which  evaluates  the  difference 
between  the  trial  FEM  model  and  the  true  damaged  FEM  model.  The  pruning  rate  is  set  to  50% 
in  this  simulation  (i.e.,  only  the  lower  50%  of  the  hypotheses  tested  defined  by  the  smallest 
errors  arc  retained  in  the  next  branching  proeess).  In  the  damage  maps,  the  elements  filled  with 
darker  eolor  have  the  higher  probabilities  of  being  damaged.  These  damage  maps  arc  created  by 
counting  the  number  of  instances  an  element  is  included  in  a  retained  trail  FEM  model;  in  effeet, 
it  is  a  histogram  of  frequency  of  inelusion  of  an  element  in  the  hypothesized  damage  set. 
Probabilistically,  the  more  frequently  an  element  appears  in  the  Bayesian  hypothesized  states,  the 
more  likely  it  is  a  truly  damaged  clement  (i.e.,  an  element  through  whieh  the  fatigue  eraek  has 
propagated).  As  shown  in  Figure  5,  the  procedure  progressively  converges  to  the  eraek  region 
and  is  terminated  when  the  variation  between  sueeessive  branching  steps  is  negligible;  for  this 
study,  this  was  found  between  the  8lh  and  911'  branching  steps. 

The  probabilistic  model-based  updating  procedure  is  also  applied  for  the  detections  of  eraek 
paths  2  and  3  as  shown  in  Figure  3(b).  These  two  eraeks  are  slightly  shorter  than  eraek  1  and 
they  are  loeated  internally  inside  the  plate.  The  detection  of  these  two  eraeks  is  more  challenging 
beeause  of  the  relatively  small  changes  in  the  modal  properties  between  the  undamaged  and 
damaged  plates  as  shown  in  Table  1.  As  shown  in  Figure  6(a),  the  model-updating  algorithm  is 
able  to  identify  aeeuratcly  the  location  of  eraek  path  2.  It  ean  be  seen  from  the  damage  map  that 
the  seareh  for  probable  damaged  elements  converges  quickly  with  a  modest  number  of  branching 
steps. 

Craek  path  3  is  a  particularly  challenging  ease  beeause  the  eraek  is  oriented  in  the  longitudinal 
direction  (i.e.,  along  the  longer  side  of  the  plate).  It  may  also  be  interesting  to  point  out  that  the 
finite  elements  are  uniformly  and  regularly  distributed  in  the  mesh.  As  shown  in  Table  1,  there 
are  very  little  ehanges  in  the  modal  properties  between  the  damaged  and  the  undamaged  plates. 
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Figure  5.  Histogram-like  damage  map  showing  the  basic  process  for  locating  probable  damage 

region  for  crack  path  1. 
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(b)  Damage  detection  for  crack  path  3 

Figure  6.  Damage  detection  results  for  individual  crack  paths  2  and  3 

Nonetheless,  as  shown  in  Figure  6(b),  the  model-updating  procedure  is  able  to  locate  the 
probable  damage  area  that  includes  the  eraek  but  the  area  is  distributed  in  the  transverse  direction 
(i.e.  along  the  shorter  side  of  the  plate).  Furthermore,  the  probable  damage  area  converges  very 
quiekly  with  very  few  branching  processes. 

Damage  Detection  of  lnJiviJual  Cracks  with  Multiple  Stiffeners:  The  model-updating  procedure 
is  applied  to  deteet  damage  in  a  plate  with  two  welded  stiffeners  as  shown  in  Figure  7(a). 
Similar  to  the  previous  eases,  the  three  eraek  paths  as  shown  in  Figure  3(b)  are  considered 
individually.  However,  the  initial  probable  damage  zone  is  expanded  to  include  the  w’eld  toes 
(i.e.,  HAZ  effeet)  around  the  two  stiffeners  and  along  the  longitudinal  edges  (i.e.,  notch  defects) 
of  the  plate.  Figure  7(b)  shows  the  probable  damage  areas  identified.  For  the  damage  detection 
for  eraek  path  1,  areas  around  the  longer  stiffener  and  along  the  top  edge  are  suspected  for 
possible  damage  very  early  in  the  process.  The  hypothetical  branches  grow  both  from  the  weld 
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(a)  A  baseline  plate  structure  with  two  stiffeners 
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(b)  Damage  detection  results  with  individual  cracks  and  multiple  stiffeners 

Figure  7.  Damage  detection  of  cracks  on  a  plate  w  ith  multiple  w  elded  stiffeners 

toe  and  edges  but  quickly  coalesce  (by  branch  step  6)  around  the  most  probable  damage  area  in 
the  vicinity  of  crack  1 .  By  branch  step  8,  the  model  updating  process  converges  and  stopped. 

For  the  detection  of  crack  2,  the  model-updating  algorithm  is  not  able  to  pin-point  the  damage 
area  with  high  probability.  Instead,  the  detection  result  suggests  two  damage  prone  areas  as 
shown  in  Figure  7(b).  Nevertheless,  the  area  around  crack  2  is  included  in  one  of  the  probable 
damage  areas.  For  detecting  the  damage  in  crack  3,  the  model-updating  procedure  identifies  the 
damage  area  successfully  with  high  probability,  as  noted  with  darker  color  elements  shown  in 
Figure  7(b).  All  in  all,  the  damage  detection  results  show  the  ability  of  the  probabilistic  model- 
updating  algorithm  in  detecting  cracks  on  a  plate  with  multiple  stiffeners  and  the  rapid 
convergence  in  identifying  the  damage  areas  even  when  the  initial  probable  damage  areas 
include  relatively  wide  spread  regions,  as  in  this  case,  around  the  stiffeners  and  along  the  edges 
of  the  plate. 

Damage  Detection  of  Multiple  Cracks:  In  this  numerical  simulation  study,  the  model-updating 
procedure  is  applied  to  detect  multiple  cracks  on  the  plate  structure.  Specifically,  two  spatially 
isolated  cracks  situated  in  the  same  plate  are  simulated  for  the  one-  (Figure  3)  and  two-stiffener 
(Figure  7(a))  hull  assemblies.  Figures  8(a)  and  8(b)  show  two  different  multi-crack  cases 
introduced  within  the  single-stiffener  plate  system.  As  shown  in  Figure  8(a),  the  model-updating 
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Multi-Crack  2 
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(a)  Crack  pattern  1,  1-stiffener 

Multi-Crack  3  Branch  5 


(b)  Crack  pattern  2,  1-stiffener 

Multi-Crack  4  Branch  6 


(c)  Crack  pattern  3,  2-stiffeners 


(d)  Crack  pattern  4,  2-stiffncrs 


Figure  8.  Damage  detection  of  multiple  cracks 


procedure  is  able  to  successfully  identify  both  cracks  next  to  the  stiffener  which  arc  located 
internally  within  the  plate.  However,  as  shown  in  Figure  8(b),  only  the  longer  and  more 
“dominant’'  crack  situated  diagonally  to  the  top  edge  of  the  plate  is  identified  by  the  damage 
detection  algorithm.  Similar  results  for  the  damage  detection  of  two  cracks  with  the  two-stiffner 
plate  system  are  shown  in  Figures  8(c)  and  8(d).  As  shown  in  Figure  8(c),  the  model-updating 
procedure  is  able  to  detect  both  internal  cracks  with  similar  size  located  next  to  each  of  the 
stiffeners.  However,  as  shown  in  Figure  8(d),  only  the  damage  from  the  longer  crack  originating 
from  the  bottom  edge  of  the  plate  is  identified.  It  is  suspected  that  the  positioning  of  a  crack  at 
the  edge  of  the  plate  may  “dominate”  the  plate  dynamics  such  that  any  craek  away  from  the  edge 
is  hard  to  identify  due  to  its  low  participation  in  the  global  plate  dynamics.  As  a  result,  further 
investigation  of  the  Bayesian  finite  element  model  updating  framework  for  the  undamaged  and 
damaged  plates  may  be  necessary  to  resolve  this  issue. 

3.2  Experimental  Simulation 

A  number  of  singlc-stiffcncr  plate  specimens  have  been  fabricated  in  full-scale  according  to  the 
schematic  shown  in  Figure  3(a).  The  intent  of  the  experimental  specimens  is  to  test  the 
applicability  of  the  probabilistic  model-updating  procedure  for  crack  damage  detection.  A  more 
detailed  overview’  of  the  experimental  setup  is  shown  in  Figure  9.  The  specimens  are  made  of 
marine  grade  aluminum  alloy  (5086).  As  shown  in  Figure  9(a),  an  18  in  long  stiffener  plate  has 
been  welded  to  each  of  the  tested  plates  with  0.625  in  long  discrete  TIG  welds  at  5  locations 
(with  a  spacing  of  4.5  in)  to  avoid  exeessive  distortions  that  might  result  from  weld  heat.  The 
base  plates  have  been  pre-hcatcd  using  a  gas  torch  to  expedite  the  welding  process  and  to 
minimize  local  warping.  As  shown  in  Figure  9(b),  the  base  plates  are  rigidly  fixed  along  the 
shorter  edges  using  8-3/8  in  stainless  screws  with  round  aluminum  washers;  furthermore,  the 
assemblage  is  rigidly  connected  to  a  large  steel  beam  post-tensioned  to  the  concrete  strong  floor 
in  the  laboratory.  Figure  9(c)  shows  a  diagonally  cut  (resembling  craek  path  2  as  shown  in 
Figure  3(b))  right  next  to  the  stiffener.  As  shown  in  the  figure,  MEMS-bascd  acecleromctcrs 
(Crossbow  CXL02LF1Z)  are  used  to  measure  the  out-of-plane  vibrations  of  the  specimen  during 
forced  vibration  testing.  Furthermore,  the  Narada  wireless  system  developed  at  the  University  of 
Michigan,  which  has  been  successfully  implemented  on  a  naval  vessel  for  hull  monitoring  ( e.g ., 
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(a)  Schematics  of  the  experimental  plate  specimens 


(b)  Plate  specimens  (c)  Specimen  with  eraek  and  accelerometers 

Figure  9.  Setup  for  experimental  tests  on  plate  specimens  with  stiffener  and  crack 


all-aluminum  FSF-1  SeaFighter),  is  deployed  as  the  primary  data  acquisition  system  for  the 
collection  of  plate  acceleration  data  (Swartz  and  Lynch  2009). 

Modal  Property  Extraction  and_  Calibration:  Both  impact  hammer  tests  (i.e.,  pure  impulsive 
loading)  and  hand  tapping  (i.e.,  colored  broadband  noise)  tests  were  applied  on  each  plate 
specimen.  For  each  plate  specimen,  albeit  the  undamaged  plate  or  the  plates  with  eraeks 
introduced,  the  hammer  test  is  applied  6  times  at  the  locations  away  from  the  weld  and  crack 
damage  area  as  shown  in  Figure  9(a).  Additionally,  each  plate  specimen  is  hand  tapped  5  times 
at  random  locations  to  simulate  colored  ambient  excitation.  The  frequency  domain 
decomposition  (FDD)  technique  is  employed  to  extract  the  modal  properties  of  the  plate 
specimens  using  the  measured  acceleration  data.  The  FDD  technique  is  widely  used  for  modal 
parameter  estimation  and  is  based  on  the  classical  complex  mode  identification  function  (Brinkcr 
et  al.  2001;  Shih  et  al.  1998).  As  illustrated  in  Figure  10,  the  estimated  modal  properties  are 
very  close  for  both  the  impact  hammer  and  hand  tapping  tests. 
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(a)  Modal  properties  extracted  from  impact  hammer  test 


Frequency  [Hz] 

(b)  Modal  properties  extracted  from  hand  tapping  test 
Figure.  10.  Estimated  modal  properties  for  the  plate  specimen  with  crack  path  1 


For  the  calibration  of  the  FEM  model,  the  modal  properties  of  the  undamaged  plate  specimen 
(without  eraeks)  arc  extracted  from  the  vibration  tests.  Slight  differences  are  observed  on  the 
modal  frequencies  between  the  experimental  results  and  the  original  finite  element  model  (used 
in  the  numerical  simulation).  Sueh  differences  arc  expeeted  beeause  of  the  variation  in  the 
idealized  finite  element  model,  plate  material  properties,  and  boundary  conditions  encountered  in 
the  experimental  setup.  To  better  calibrate  between  the  numerical  model  and  the  experimental 
specimens,  the  elastie  (Young’s)  modulus  used  in  the  finite  clement  model  is  redueed  to  90%  of 
the  nominal  values  (10,300  ksi)  so  that  the  first  mode  frequencies  of  the  undamaged  plate  from 
both  the  FEM  model  and  the  experimental  tests  are  matehed.  Table  2  summarizes  the  modal 
frequencies  of  the  plate  structure  with  eraek  path  1  obtained  from  the  experimental  hammer  and 
hand  tapping  tests  as  well  as  from  the  updated  finite  element  model.  It  ean  be  seen  that  the 
results  mateh  with  the  modal  assuranee  criteria  (MAC)  exceeding  90%. 

Damage  Identification  Results ;  The  model  updating  procedure  is  applied  to  identify  the  damage 
area  (eraeks)  based  on  the  modal  properties  extracted  from  the  plate  specimens.  Figure  1 1  shows 
the  damage  identification  results  for  the  three  eraek  paths  from  both  impaet  hammer  and  hand 
tapping  tests.  It  ean  be  seen  that  the  model  updating  procedure  successfully  deteets  eraek  paths  1 
and  3.  However,  for  eraek  path  2,  although  the  most  probable  damage  area  is  identified  near  the 
damage  area,  the  procedure  fails  to  pin-point  the  aetual  eraek  location.  It  ean  also  be  observed 
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Table  2.  Summary'  of  modal  properties  extracted  for  specimen  with  eraek  path  1 


Impact  hammer  test 

Modal  frequency 

Mode  shape 

Mode 

Test 

Updated  FEM 

MAC 

(Hz) 

(Hz) 

1 

24.22 

25.17 

0.998 

2 

35.89 

37.89 

0.983 

3 

65.43 

67.36 

0.992 

4 

85.57 

85.14 

0.982 

5 

122.07 

120.43 

0.932 

Hand-tapping  test 

Modal  frequency  Mode  shape 


Mode 

Test 

Updated  FEM 

MAC 

(Hz) 

(Hz) 

1 

24.22 

25.17 

0.999 

2 

35.94 

37.89 

0.996 

3 

65.43 

67.36 

0.985 

4 

85.64 

85.14 

0.976 

5 

121.97 

120.43 

0.935 

Target  Impact  hammer  test  Hand  tapping  test 
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Figure  11.  Damage  identification  results  from  experimental  test  specimens 


that,  in  general,  the  damage  identification  results  based  on  impact  hammer  tests  are  more  precise, 
as  illustrated  from  the  darker  colors  shown  in  the  damage  maps.  It  should  be  cautioned  that,  as 
in  any  experimental  test,  the  fabrication  qualities  of  the  plate  specimens  and  the  welds  do  vary. 
The  modal  frequencies  could  vary  up  to  5-7%.  Further  study  on  the  impact  of  modeling 
uncertainties  is  needed.  The  effects  of  noise  on  the  probabilistic  model-updating  procedure  are 
also  worth  further  investigation. 
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4.  Summary  and  Future  Work 

A  model-based  damage  detection  algorithm  has  been  explored  for  structural  health  monitoring  of 
aluminum  ship  hulls.  A  particular  emphasis  was  placed  on  the  detection  of  fatigue-related  crack 
damage  in  plate-type  hull  assemblies.  The  algorithm  combines  vibration-based  damage 
detection  with  a  finite  element  method  model  to  provide  the  basis  for  high-fidelity  damage 
detection  using  a  limited  number  of  sensors  installed  in  the  assembly.  The  FEM  model  of  the 
interested  structure  is  successively  updated  in  a  Bayesian  inference  scheme  employing  a  Braneh- 
and-Bound  technique  as  the  primary  systematic  searching  strategy.  Numerical  simulations  and 
experimental  tests  have  been  conducted  to  illustrate  the  model-based  monitoring  framework  and 
to  validate  the  potential  application  of  the  probabilistie-based  damage  detection  approach. 
Preliminary  numerical  and  experimental  results  indicate  that  the  model-based  damage  detection 
procedure  is  able  to  successfully  identify  crack  damage  areas,  especially  when  the  size  of  the 
crack  is  relatively  large  and  strongly  influences  the  global  plate  dynamics. 

Currently,  the  Bayesian  model  updating  framework  is  being  applied  to  a  more  complex,  but 
realistic  hull  component.  Specifically,  a  scries  of  fatigue  tests  of  the  MAHI  (Monitored 
Aluminum  Hull  Integrity)  specimens  (welded  aluminum  assemblages  of  a  bulkhead  element  and 
wide  flange  column)  are  underway  at  the  Naval  Surface  Warfare  Center-Cardcroek  Division 
(NSWCCD).  Our  research  efforts  have  included  active  engagement  with  researchers  at 
NSWCCD  on  the  fatigue  tests  both  in  instrumentation,  data  collection  and  finite  clement 
modeling  of  the  specimens.  Our  future  plan  is  to  further  apply  and  validate  the  model-based 
damage  detection  framework  using  the  MAHI  project. 

Frequent  inspection  of  ship  hulls  can  extend  the  operational  life  of  a  ship  since  the  detection  of 
the  onset  of  damage  can  reduce  overall  ship  lifc-eyele  costs.  However,  current  visual  inspections 
of  the  entire  hull  are  both  costly  and  labor-intensive.  Therefore,  the  instrumentation  of  a 
structural  health  monitoring  (SHM)  system  coupled  with  an  effective  damage  detection 
methodology  such  as  the  approach  presented  in  this  report  can  reduce  the  cost  of  inspection  by 
providing  inspectors  with  a  prioritized  list  of  probable  areas  of  damage.  Such  a  list  would 
effectively  narrow  down  the  areas  of  the  hull  requiring  detailed  inspections  (perhaps  by 
conducting  nondestructive  testing  with  ultrasonic  waves,  etc.)  and  repairs.  Our  future 
investigation  will  extend  the  current  damage  detection  methodology  as  an  integral  component  of 
a  life  cycle  design,  monitoring  and  maintenance  framework  for  ship  structure. 
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